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Abstract 

A real photon has a complicated nature, whereby it may remain unresolved 
or fluctuate into a vector meson or a perturbative qq pair. In 77 events, this 
gives three by three combinations of the nature of the two incoming photons, 
and thus six distinct event classes. The properties of these classes are partly 
constrained by the choices already made in our related 7p model. It is 
therefore possible to predict the energy- dependence of the cross section for 
each of the six components separately. The total cross section gives support 
to the idea that a simple factorized ansatz with a pomeron and a reggeon 
term can be a good approximation. Event properties undergo a stepwise 
evolution from pp to 7p to 77 events, with larger charged multiplicity, more 
transverse energy flow and a higher jet rate in the latter process. 
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1 Introduction 



There are many reasons for being interested in 77 physics. The coUision between two 
photons provides the richest spectrum of (leading-order) processes that is available for 
any choice of two incoming elementary particles. For instance, since the photon has a 
hadronic component, all of hadronic physics is contained as a subset of the possibilities. 
Additionally, the photon can appear as an unresolved particle or as a perturbative qq fluc- 
tuation, giving a host of possible further interaction processes. The relative amount of 
these components and their respective properties are not fully understood today. A correct 
description of the components of the total 77 cross section and the related event shapes 
therefore is the ultimate challenge of 'minimum-bias' physics. Specific issues include the 
description of the photon wave function, duality between perturbative and nonperturba- 
tive descriptions of the resolved photon, the role of multiple parton-parton interactions 
and the related minijet phenomenology and eikonalization of the total cross section, the 
transition between soft and hard physics, the transition between the real photon and the 
virtual one, and so on. 

In addition to the direct reasons, there are also indirect ones. The process e^e~ 
e"'"e~77 e"'"e~X will be a main one at LEP 2 and future linear e^e~ colliders. Therefore, 
77 events are always going to give a non-negligible background to whatever other physics 
one is interested in. The devising of efficient analysis strategies must be based on a good 
understanding of 77 physics. 

The study of 77 physics has a long history, and it is not our intention here to give 
a complete list of references. Many topics have been covered by contributions to past 
workshops In recent years HERA has been providing rich information on the related 
7P processes, and thereby stimulating the whole field 0. Further developments can be 
expected here. Currently 77 interest is focussed on LEP 2 P|. Already LEP 1.5 has 
provided ample reminder of the important role of 77 processes when away from the Z° 
pole. LEP 2 brings the promise of a large event rate at reasonably large 77 energies. In 
the future, the laser backscattering option of linear e+e~ colliders offers the promise of 
obtaining a large rate of very high-energy 77 interactions, typically at up to 70% of the 
energy of the corresponding e^e~ collisions. Then the two aspects above come together 
in force, both with new chances to understand photon interactions and new challenges to 
eliminate the 77 background to other processes. 

The starting point for the current paper is our model for 7p physics [Q. Many of 
the basic assumptions can be taken over in an (almost) minimal fashion, while further 
new ones appear. Based on the experience from HERA, some old assumptions can be 
sharpened and further developed. LEP 2 and future linear colliders will allow new tests 
to be carried out. Our recent studies on the parton distributions of the photon are 
parts of the same physics program, and provide further important building blocks for the 
current study. In this paper we also emphasize the gradual evolution from pp to 7p to 77 
events, that allows some cross-checks to be carried out systematically. Parts of this work 
has already been presented in a preliminary form at workshops P, |^. 

No model exists in a vacuum. For the approach we are going to take, one important 
line of work is the subdivision of photon interactions by the nature of the photon P]. 
Minijet phenomenology has attracted much attention in recent years Other related 
works will appear as we go along. However, none of these approaches attempts to give a 
complete description of 77 cross sections and event properties, but only concentrate on 
specific topics. Here we will try to be more ambitious, and really provide all the necessary 



1 



aspects in one single framework. The only other work with a somewhat similar global 



scope is the recent studies within the context of dual topological unitarization . 

Some main areas are still left out of our description. In all that follows, both incoming 
photons are assumed to be on the mass shell; further issues need to be addressed when 
either photon or both of them are virtual. The issue of the eikonalization of the anomalous 
and direct components of the photon wave function will be partly deferred to future 
studies; the evidence for some form of eikonalization will become apparent as we go 
along. Finally, for reasons of clarity, we restrict ourselves to discussing what happens in 
the collision between two photons of given momenta. The addition of photon flux factors 
0] complicates the picture, but does not add anything fundamentally new. 

Section 2 contains a description of the photon wave function and 77 event classes, 
section 3 of total and partial cross sections, section 4 of event properties and section 5 
gives a summary and outlook. 



2 The Photon Wave Function and Event Classes 

To first approximation, the photon is a point-like particle. However, quantum mechani- 
cally, it may fiuctuate into a (charged) fermion-antifermion pair. The fiuctuations 7 qq 
are of special interest to us, since such fiuctuations can interact strongly and therefore 
turn out to be responsible for the major part of the 7p and 77 total cross sections, as we 
shall see. On the other hand, the fiuctuations into a lepton pair are uninteresting, since 
such states do not undergo strong interactions to leading order, and therefore contribute 
neghgibly to total hadronic cross sections. The leptonic fiuctuations are perturbatively 
calculable, with an infrared cut-off provided by the lepton mass itself. Not so for quark 
pairs, where low-virtuality fiuctuations enter a domain of non-perturbative QCD physics. 
It is therefore customary to split the spectrum of fiuctuations into a low-virtuality and 
a high-virtuality part. The former part can be approximated by a sum over low- mass 
vector-meson states, customarily (but not necessarily) restricted to the lowest-lying vec- 
tor multiplet. Phenomenologically, this Vector Meson Dominance (VMD) ansatz turns 
out to be very successful in describing a host of data. The high-virtuality part, on the 
other hand, should be in a perturbatively calculable domain. 
In total, the photon wave function can then be written as 

I7) = Cbare|7bare) + E + ^ Cq|qq) + H (1) 

V=p'^ ,LO,<f>,i q=u,d,s,c,b £=e,/i,T 

(neglecting the small contribution from T). In general, the coefficients Cj depend on the 
scale /i used to probe the photon. Thus c\ ~ (aom/2vr)(2/3) ln(/x^/m^). Introducing a cut- 
off parameter fco to separate the low- and high-virtuality parts of the qq fiuctuations, one 
similarly obtains Cq ~ (aem/27r)2eqln(/i^//co). Since each qq fiuctuation is characterized 
by some virtuality or transverse momentum scale /c, the notation in eq. (|1]) should really 
be viewed as shorthand, where the full expression is obtained by 

ej«)_^2</;'|!|«;.^), (2) 

and correspondingly for the lepton component. This is the form assumed in the following. 
The VMD part corresponds to the range of qq fiuctuations below /cq and is thus /i- 
independent (assuming > ko). In conventional notation Cy = inaem/fvy with /y/47r 
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determined from data to be 2.20 for p°, 23.6 for uj, 18.4 for cj) and 11.5 for J/-?/; |TT|. 
Finally, Cbare is given by unitarity: c^^^^ = ^3 = 1 - E Cy - E Cq - E c|. In practice, 
Cbare IS always close to unity. Usually the probing scale is taken to be the transverse 
momentum of a 2 — > 2 parton-level process. Our fitted value ~ 0.5 GeV @] then sets 
the minimum transverse momentum of a perturbative branching 7 qq. 

In part, fco is an unphysical parameter, and one would expect a continuity under 
reasonably variations of it. That is, if the VMD sum were to be extended beyond the 
lowest-lying vector mesons to also include higher resonances, it should be possible to com- 
pensate this by using a correspondingly higher cut-off for the continuous perturbative 
spectrum. This may provide some guidelines when exploring the physics implications of 
the ansatz above. The VMD-perturbative state duality has its limits, however. Higher 
excitations of vector mesons have larger wave-function radii than the lowest-lying states 
when each is produced 'on shell' in the time-like region (roughly r oc m), while the uncer- 
tainty relation gives smaller radii for the higher-virtuality components of the real photon 
wave function (r oc 1/m). Correspondingly, the contributions of u and could not 
well be described by perturbation theory alone: the cy parameters are related to the 
absolute rates of the elastic processes 7p Vp, and here the observed relation between 
p'' and uj production is in agreement with the 9 : 1 VMD expectations of coherent vector 
meson wave-functions, while incoherent interactions of perturbative components |uu) and 
|dd) would have lead to equal production of p° and u. 

The subdivision of the above photon wave function corresponds to the existence of 
three main event classes in 7p events, cf. Fig. |l]: 

1. The VMD processes, where the photon turns into a vector meson before the inter- 
action, and therefore all processes allowed in hadronic physics may occur. This in- 
cludes elastic and diffractive scattering as well as 1ow-j9_l and high-p_L non-diffractive 
events. 

2. The direct processes, where a bare photon interacts with a parton from the proton. 

3. The anomalous processes, where the photon perturbatively branches into a qq pair, 
and one of these (or a daughter parton thereof) interacts with a parton from the 
proton. 

All three processes are of 0(Q;em)- However, in the direct contribution the 'parton' distri- 
bution of the photon is of 0(1) and the hard scattering matrix elements of 0{aem)i while 
the opposite holds for the VMD and the anomalous processes. As we already noted, the 
fluctuations are not interesting, and there is thus no class associated with them. 

The above subdivision is not unique, or even the conventional one. More common is 
to lump the jet production processes of VMD and anomalous into a class called resolved 
photons. The remaining 'soft- VMD' class is then defined as not having any jet production 
at all, but only consisting of low-p^ events. We find such a subdivision counterproductive, 
since it is then not possible to think of the VMD class as being a scaled-down version (by 
a factor Cy) of ordinary hadronic processes — remember that normal hadronic collisions 
do contain jets part of the time. 

In a complete framework, there would be no sharp borders between the three above 
classes, but rather fairly smooth transition regions that interpolate between the extreme 
behaviours. However, at our current level of understanding, we do not know how to 
do this, and therefore push our ignorance into parameters such as the k^ scale and the 
fy/^T^ couplings. From a practical point of view, the sharp borders on the parton level 
are smeared out by parton showers and hadronization. Any Monte Carlo event sample 
intended to catch a border region would actually consist of a mixture of the three ex- 
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treme scenarios, and therefore indeed be intermediate. This issue is discussed further in 



section |373 . 

The difference between the three classes is easily seen in terms of the beam jet struc- 
ture. The incoming proton always gives a beam jet containing the partons of the proton 
that did not interact. On the photon side, the direct processes do not give a beam jet 
at all, since all the energy of the photon is involved in the hard interaction. The VMD 
ones (leaving aside the elastic and diffractive subprocesses for the moment) give a beam 
remnant just like the proton, with a 'primordial k±^ smearing of typically up to half a 
GeV. The anomalous processes give a beam remnant produced by the 7 — >■ qq branching, 
with a transverse momentum going from ko upwards. Thus the transition from VMD to 
anomalous should be rather smooth. 

A generalization of the above picture to 77 events is obtained by noting that each of 
the two incoming photons is described by a wave function of the type given in eq. (|lD. In 
total, there are therefore three times three event classes. By symmetry, the 'off-diagonal' 
combinations appear pairwise, so the number of distinct classes is only six. These are, cf. 
Fig. |: 

1. VMD X VMD: both photons turn into hadrons, and the processes are therefore the 
same as allowed in hadron-hadron collisions. 

2. VMD X direct: a bare photon interacts with the partons of the VMD photon. 

3. VMD X anomalous: the anomalous photon perturbatively branches into a qq pair, 
and one of these (or a daughter parton thereof) interacts with a parton from the 
VMD photon. 

4. DirectX direct: the two photons directly give a quark pair, 77 — > qq. Also lepton 
pair production is allowed, 77 , but will not be considered by us. 

5. DirectX anomalous: the anomalous photon perturbatively branches into a qq pair, 
and one of these (or a daughter parton thereof) directly interacts with the other 
photon. 

6. Anomalous X anomalous: both photons perturbatively branch into qq pairs, and 
subsequently one parton from each photon undergoes a hard interaction. 

The first three classes above are pretty much the same as the three classes allowed in 7p 
events, since the interactions of a VMD photon and those of a proton are about the same. 
The main parton-level processes that occur in the six classes are: 

• The 'direct' processes 77 qq only occur in class 4. 

• The '1-resolved' processes 7q qg and 7g qq occur in classes 2 and 5. 

• The '2-resolved' processes qq' qq' (where q' may also represent an antiquark), 
qq ~^ q'q') qq ~^ gg^ qg ~^ qg) gg ~^ qq and gg — > gg occur in classes l, 3 and 6. 

• Elastic, diffractive and low-p± events occur in class 1. 

In the list above we have only indicated the lowest-order processes. Within the context 
of the leading-log approximation, at least, the subdivision into six event classes is easily 
generalized to graphs with an arbitrary number of partons in the final state. This clas- 
sification is illustrated in Fig. |^ for a generic ladder graph. To this picture final-state 
radiation can be added trivially. 

The notation direct, 1-resolved and 2-resolved is the conventional subdivision of 77 
interactions. The rest is then called 'soft- VMD'. As for the 7p case, our subdivision is 
an attempt to be more precise and internally consistent than the conventional classes 
allow. One aspect is that we really want to have a VMD x VMD class that is nothing 
but a scaled-down copy of the p°p° and other vector-meson processes, with a consistent 
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transition between low-p± and high-p^ events (see below). Another aspect is that, in a 
complete description, the VMD and anomalous parts of the photon give rise to different 
beam remnant structures, as discussed above, even when the hard subprocess itself may 
be the same. 

A third aspect is that our subdivision provides further constraints; these, at least in 
principle, make the model more predictive. In particular, the parton distributions of the 
photon are constrained by the ansatz in eq. (|^) to be given by 

Paix, /i') = /.^) + /i^) + /x^; kl) . (3) 

Here 
and 

The anomalous part, finally, is fully calculable perturbatively, given the boundary condi- 
tion that the distributions should vanish for /x^ = k^. In principle, everything is therefore 
given. In practice, the vector-meson distributions are not known, and so one is obliged 
to pick some reasonable ansatz with parameters fitted to the data. This is the approach 
taken in our SaS parameterizations 0. By comparison, conventional distributions are 
defined for resolved processes only: 

f2''\^, = /i^) + /i^; ki) . (6) 

These resolved distributions are then less constrained, in particular with respect to the 
momentum sum of resolved partons. 



3 Cross Sections 

3.1 The total cross section and its subdivision 

Total hadronic cross sections show a characteristic fall-off at low energies and a slow rise 
at higher energies. This behaviour can be parameterized by the form 

at^{s)=X^^s^ + Y^^s-^ (7) 



for A ^ B ^ X . The powers e and t] are universal, with fit values 

e ^ 0.0808 , T] « 0.4525 , (8) 

while the coefficients X^^ and Y^^ are process-dependent. Equation (|^ can be in- 
terpreted within Regge theory, where the first term corresponds to pomeron exchange 
and gives the asymptotic rise of the cross section. Ultimately, this increase violates the 
Froissart-Martin bound W^:, e should therefore be thought of as slowly decreasing with 



energy (owing to multi-pomeron exchange effects), although data at current energies are 
well fitted by a constant e. The second term, the reggeon one, is mainly of interest at 
low energies. For the purpose of our study we do not rely on the Regge interpretation of 
eq. (1^, but can merely consider it as a convenient parameterization. 
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The VMD part of the 7p cross section should have a similar behaviour. The direct 
part reflects the parton distributions of the proton; a small-x behaviour like xf{x) ~ 
would give cTj? ~ s^ The anomalous part is less easily classified: a purely perturbative 
description would not give a behaviour like VMD, but a duality argument with anomalous 
states interpreted in terms of higher vector-meson states would. Empirically, the 7p data 
are well described by 

aZis)^Q7.7s' + 129s-^ [/ib], (9) 

with s in GeV^. (Cross-sections are throughout given in mb for hadron-hadron interac- 
tions, in /ib for 7-hadron ones and in nb for 77 ones.) Actually, the above formula is 
a prediction |13| preceding the HERA data [15, The conclusion would seem to be 



that, at least as far as total cross sections are concerned, the extended VMD description 
of anomalous interactions is a reasonable first approximation. 

If we then take the Regge-theory ansatz seriously also for the photon, it is possible to 
derive an expression for the total 77 cross section 

~ 2115^ + 215 5-" [nb]. (10) 

This is based on the assumption that the pomeron and reggeon terms factorize, X"^^ = 
PajpPbjp and Y^^ = 7air7bir, so that X^^ = {X^^f/X^P and Y^'' = 2{Y^Pf /{Ypp+Y^p), 
with ^ 21.70 and {Y^^ + Y^^)/2 ^ 77.23. In hadronic cross sections, factorization 
seems valid for the pomeron term but not for the reggeon one, e.g. X^p = Xpp while 
ypp PS 98.39 > ypp ^ 56.08. The choice of using the average of Y^p and Fpp then is an 
arbitrary one, though it can be motivated roughly by arguments of counting the number of 
allowed valence quark/ ant iquark annihilation/exchange diagrams possible in the various 
processes. The band of uncertainty can be obtained by using either ypp or Y^^ alone, 
i.e. Y'^'^ = 297 and 169. This ambiguity only affects the low-energy behaviour, and so is 
not critical for us. It is illustrated in Fig. El where we also compare with existing data on 



Note that factorization is assumed to hold separately for the pomeron and the reggeon 
terms, not for the total cross section itself. That is, the relation = 2((Tto't)^/(o"tot + 
atot) is not exact in this approach, although numerically it is a very good approximation 
(usually to better than 1%). 

Our eq. (0) above should be compared with the time-honoured expression a"''^ = 
240 + 270/W [0. This corresponds to a critical pomeron, e = 0, as was commonly 
assumed in the early seventies, and an 77 = 0.5, but it is otherwise in the same spirit as 
our formula. Also numerically the two closely agree at not too large energies, see Fig. ^ 

One should remember that our expression ( p!OD is here 'derived' based on a simple 
Regge-theory ansatz that has no real validity for the photon. Next we will proceed to 
study the contributions of the individual event classes. The constraints that come from 
7P physics data then directly feed into constraints on the contribution from these classes 
and therefore on the total 77 cross section. At the end of the day we will therefore show 



that a cross section behaving roughly like eq. (|T0D should be a good approximation. In 
doing so, the properties of the event classes are also fixed, to a large extent. 

Based on the subdivision into event classes, the total 7p cross section may be written 

as 

and the total 77 one as 

Otot — OvMDxVMD ^ ^"VMDxdir ^ ^'^VMDxanom ^ '^dirxdir ^ ^'^dirxanom ^ Oanomxanom • \^'^) 
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Here we explicitly keep the factor of 2 for the off-diagonal terms, where the role of the 
two incoming photons may be interchanged. 



3.2 The VMD contributions 

The Vp cross sections may be parameterized as 

af^r ^ ^ro"? ^ I {ctC'' + afo?) ^ 13.63 + 31.79 s'^ [mb], (13) 
crt ^ ^fl" + o^fo/ - <^ ^ 10.01 - 1.51 s-" [mb]. (14) 
The 0p cross section is not expected to have a reggeon term and indeed the additive 



quark model formulae give a contribution close to zero; a small negative term could 
easily come from threshold effects and so we choose to keep it. Lacking measurements 
of Dp cross sections we cannot use the additive quark model to estimate the J/V^p cross 
section. The latter could in principle be extracted from data on^/ip production in nuclear 
collisions. Superficially a value of about 6 mb at 5 10 GeV comes out but this value 

presumably is too large since in the so far accessible kinematical range of large positive 
xp the J/?/' is formed outside the nucleus. Therefore we fix the VMD coupling of the J/V' 
at its leptonic value and use a low-energy measurement of elastic ^/ip photoproduction 
to determine the J/^/'p cross section. This implies a reduction of the V'P cross section by 
a factor of about 10 compared to the 0p one, in agreement with the expectation that the 
soft pomeron couples more weakly to heavier quarks. Again using factorization for the 
pomeron and reggeon terms separately, the total cross section for two vector mesons is 

^tot ' ~ + ^ s . (15) 

These X and Y coefficients are collected in Table |1|. 

The total VMD cross sections are obtained as weighted sums of the allowed vector- 
meson states, 

<MD = E^^^t''o?^54.^ + 115s-'' [/ib], (16) 

V Jv 

7mdxvmd = E ^ E ^ <^r^' - 133 + 170 [nb]. (17) 

VI JVl V2 JV2 



a 



In Fig. 1^ we show the breakdown of ctvmdxvmd by vector-meson combination. Obviously 
the combination dominates. 

For a description of VMD events, a further subdivision into elastic (el), diffractive 
(sd and dd for single and double diffractive) and non- diffractive (nd) events is required. 
Keeping only the simplest diffractive topologies, one may write 

4o^{s) = ai\s) + atH^^^is) + + 4i{s) + ati{s) . (18) 

The elastic and diffractive cross sections for all required Vp and V1V2 processes have 
been calculated and parameterized in the context of our model presented in ref. pO| . 
The same formulae are used as those collected in section 4 of that paper, and so are not 
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repeated here; only the expressions in its eq. (26) have to be replaced. The following 
parameterizations have been chosen: 





= CiS + C2 , 




BxB 


= C3 H , 

s 




-^^-^ max, AX 


= C5S + C6 , 




Bax 


= C7 + — , 
S 




Ao 


d2 

= di + — + 
m s 


ds 

1 2 ' 
m s 


^max,XB 


= S d4 + - — 

\ m s 


de 


Bxx 


= dj -\ — -= + 


dg 
S 



(19) 

The coefficients Cj and di are given in Table |l|. Additionally the b slope parameters are 
bp = 2.3 GeV-2, bp = b^ = b^ = 1.4 GeV-^ and fej/^ = 0.23 GeV-^. The non-diffiactive 
cross-section is then given by whatever is left. This subdivision is shown in Fig. ^ for the 
sum of all meson combinations, which then mainly reflects the p^p^ composition. 

The (Tnd may be further subdivided into a low-p_L and a high-p^ class. Since the 2 — > 2 
parton-parton scattering cross sections are divergent in the limit p± — 0, some further 
care is needed for this classification. We expect the perturbative formulae to break down 
at small since an exchanged gluon with a large transverse wavelength A_l ~ l/pi_ 
cannot resolve the individual colour charges inside a hadron. The hadron being a net 
colour singlet, the effective coupling should therefore vanish in this limit. A parameter 
P±min = P±min(s) is iutroduccd to dcscribc the border down to which the perturbative 
expression is assumed to be valid The jet rate above p_Lmin may still be large, in fact 
even larger than the total and- It is therefore necessary to allow for the possibility of 
having several perturbative parton-parton interactions in one and the same event, i.e. to 



unitarize the jet emission probability. We do this using the formalism of ref. A fit 

to collider multiplicities gives 



Px.in(.) = pl'Zis) - 1.30 + 0.15 ^[1^ [GeV] . (20) 



Here we uses the CTEQ 2L leading-order parton distributions (extended to small x 



and Q as described in [Q), with p]_ as scale choice. 



3.3 The direct and anomalous contributions 

Comparing eqs. (H) and ([T6|) , about 80% of the 7p total cross section is seen to come 
from the VMD term. The remaining 20% is to be attributed to the direct and anomalous 
components. When applying a perturbative description, the anomalous part is negligible 
at small energies. The dependence of the direct cross section on /cq can then be used to 
determine this parameter. We obtain a value of ~ 0.5 GeV Q], which is consistent 
with the simple-minded answer ~ In our study of the parton distributions of 

the photons a reasonable f2''^'^^{x, p?) was obtained with Qq = 0.6 GeV, i.e. the same 
order. For this study we have stayed with the latter number, /cq = 0.6 GeV. 
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The anomalous process contains two cut-off parameters, the fco scale for the photon to 
branch to a perturbative qq pair and a Pj^^ scale for one of the anomalous-photon partons 
to interact in a hard process. As a first guess, one might choose also to be given by 

eq. (pOl). However, this turns out to give a cross section increasing too rapidly at large 
energies. Physically, it is understandable why hard processes should be more suppressed 
at small p± in anomalous processes than in VMD ones: the anomalous photon corresponds 
to a qq pair of larger virtuality than a VMD one, and hence of smaller spatial extent, i.e. 
with larger potential for colour screening. The best recipe for including this physics aspect 
is not well understood. As a purely pragmatical recipe, one can pick Pji^^{s) with an s 
dependence such that the VMD, direct and anomalous processes add up to the expected 
behaviour (P) . Over the energy range 20 1000 a suitable parameterization then is 

PTZnis) ~ 0.6 + 0.125 ln'(l + v^/10) [GeV] . (21) 

This is based on parton distributions SaS ID for the photon and CTEQ 2L for the 
proton, combined with lowest-order matrix elements. At low energies the results are fairly 
unstable to variations in kQ, and so the behaviour of p±^^{s) in this region should not 
be over- interpreted. The split of the 7p cross section into VMD, direct and anomalous is 
shown in Fig. |^. 

In this purely perturbative picture we have neither included any class of 'soft' anoma- 
lous interactions nor the possibility of multiple parton-parton interactions. Keeping ev- 
erything else the same, the former would increase the cross section and the latter decrease 
it. However, when introducing a soft component, it is important to avoid double-counting. 
To illustrate the issue, consider the simple graph of Fig. There are two transverse 
momentum scales, kj_ and pj_. It is a simpler version of Fig. |^, with inessential gluons 
removed and for 7p rather than 77, so has one scale less. The allowed phase space can 
then conveniently be represented by a two-dimensional plane. Fig. ^b. The region kj_ < kQ 
corresponds to a small transverse momentum at the 7 — > qq vertex, and thus to VMD 
processes. For k± > ko, the events are split along the diagonal k± = p±- If kj_ > p±, the 
hard 2 — 2 process of Fig. ||a is 7g — qq, and the lower part of the graph is part of the 
leading log QCD evolution of the gluon distribution inside the proton. These events are 
direct ones. If p± > k±, on the other hand, the hard process is qq' — »• qq', and the 7 ^ qq 
vertex builds up the quark distribution inside a photon. These events are thus anomalous 
ones. 

By analogy with the VMD representation, each fixed-fc_L component of 7 <-> qq fluctua- 
tions can be considered as a separate 'hadron species', with a density of states proportional 
to dk'^/kj_. Each vertical 'tower' at some given k± scale would correspond to a higher 
excited vector resonance in the context of a generalized VMD model. In this tower, the 
soft events would be in the direct sector and the hard events in the anomalous sector. In 
the region of large k_i_ values the perturbative language is well defined, and no problems 
should arise. As smaller and smaller k±s are considered, however, one could expect event 
properties that are intermediate to those of VMD. In particular, multiple parton-parton 
interactions could be possible, and this would affect the relation between calculated jet 
cross sections and the total event cross section. Previously we had to introduce a large 
PTndn scale at high energies to solve the problem of too large an anomalous cross section, 
which means we left an un-populated hole in the middle of Fig. (indicated by a ques- 
tion mark). The hope is that multiple interactions would provide a natural resolution of 
this problem, in the sense that most anomalous events have one hard scattering above 
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PTmhi' while the anomalous region with p± < does not significantly contribute new 

events but rather further interactions inside the events above. 

The hope for a simple solution is not borne out by studies, however . Eikonalization 
does dampen the increase of the anomalous and direct cross sections, but not enough. 
With a Px'^is) = p±^ii{s) the cross section is still increasing too fast at high energies, if 
the pomeron-style behaviour is taken as reference. Problems that appear already in the 7p 
sector are even more severe in attempts at a 77 description along the same lines. It there- 
fore seems clear that further aspects have to be taken into account, such as momentum 
conservation, coherence effects or strict geometrical cuts. 

We intend to return to these problems, but for the moment stay with a purely pertur- 
bative description of the direct and anomalous components. By pushing this approach to 
its logical conclusion, we see what to expect from it and what limitations it has. Further- 
more, for most event properties we expect the perturbative description to be perfectly 
adequate. 



3.4 The total 77 cross section by component 

Turning to the 77 cross sections, in principle all free parameters have now been fixed, 
and the cross section for each of the six event classes can be obtained. The VMDxVMD 
one has already been discussed; the others are given as integrals of 2 — > 2 scattering 
cross sections above the respective p± cut-offs already specified. The results are shown in 
Fig. 1^, class by class. For comparison, we also show the results that would be obtained 
in a simple factorization ansatz 



1 _L A. . ^PP I ^PP 



CpTT. = Z "J (22) 



where i,j = VMD, direct and anomalous. 
A few comments about each of the classes: 

1. For the VMDxVMD class in principle the simple factorization ansatz is exact in 
our model; some minor deviations come from the reggeon term. 

2. The VMD X direct class comes out about a factor 3/2 larger than expected from the 
factorized ansatz. This difference can be understood by comparing the jet and the 
total cross sections of a proton and a vr meson, where the latter is taken as a proto- 
type for a VMD meson. Both the p and the vr have parton distributions normalized 
to unit momentum sum, and the same small-x behaviour (using our prescription 
[^]). Neglecting some differences in the shape of the parton distributions, the jet 
rates therefore are comparable between pp, vrp and tttt collisions. The total cross 
sections, on the other hand, are scaled down roughly by a factor 2/3 between pp 
and vrp. Therefore the jet rate per event is a factor 3/2 larger for vrp than for pp, 
and it is this factor that appears above. That is, eq. (l2^) would have worked only 
if 77r/7r7r cross sections could have been used rather than 7p/pp ones. The larger 
jet rate per event for mesons should be reflected in differences in the eikonalization 
treatment of the direct and anomalous components of 7p and events. 

3. The VMD X anomalous component gives exactly the same factor 3/2 mismatch as 
discussed above for the VMD x direct one. 

4. The direct X direct component is not at all well predicted by the factorized ansatz. 
The latter yields a cross section growing at large energies at a rate related to the 
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small-a; behaviour of the proton distribution functions, i.e. oc for our modified 
distributions. On the other hand, the total cross section for 77 — * qq is proportional 
to ln(s/fcg)/s, and thus drops rapidly with cm. energy. 

5. The direct X anomalous component compares reasonably well with the prediction 
from factorization. 

6. The anomalous X anomalous process, finally, is most uncertain, since it completely 
involves the interactions of the least well understood component of the photon wave 
function. 

In Fig. |10| the total 77 cross section is compared between the Regge type ansatz ([T0|) 
and the sum of the six classes above, eq. (0). It should be remembered that the first 
three are the dominant ones. In fact, since the direct and anomalous components to- 
gether give about 20% of the 7p total cross section, the expectation is that the last 
three classes together would only give a 4% contribution to the total 77 cross section. 
The anomalous X anomalous component may give a somewhat larger contribution than 
expected, but still the 4% number gives the right ballpark. The first three classes, on the 
other hand, are all related to the respective 7p classes, with only a replacement of a p 
by a y (and an extra weight factor 47raem//y)- Apart from the appearance of a factor 
3/2 in the VMDx direct and VMD x anomalous components, which should (largely if not 
completely) go away in a fully eikonalized description, these components behave as ex- 
pected. This makes the argumentation for eq. ([T0|) credible. However, if one wants to take 
a conservative approach, the spread between the two curves in Fig. [TO could be viewed 
as a band of uncertainty. The data are not yet precise to provide any discrimination, cf. 
Fig. |. 

One can also compare our a'^J^ with the numbers obtained in various minijet-based 
approaches [g i. For E^^ = 200 GeV, cross sections in the range 1000-1800 nb are 
typically obtained, but are reduced to about 500 nb if unitarity is enforced, in agreement 
with our results. 



4 Event Properties 

The subdivision of the total 7p and 77 cross sections above, with the related choices of 
cut-off parameters etc., specifies the event composition at the hard-scattering level. Some 



interesting observables can be based on this classification alone. For instance. Fig |TT] gives 
the cross section for elastic events of the kinds p°a;, p'^cf) and p^J/ip. The former 

three processes are good tests for the validity of the VMD ansatz, whereas the last one 
could provide new insights. 

For most studies it is necessary to consider the complete event structure, i.e. to add 
models for initial- and final-state QCD radiation (parton showers), for beam remnants, 
and for fragmentation and secondary decays [Q. A Monte Carlo generation of complete 
hadronic final states is obtained with Pythia/Jetset |2^. Thus any experimental quan- 



tity can be studied. This section gives some representative examples. In particular, we 
compare the properties of pp, 7p and 77 events. It should be noted that pp and pp 
events are very similar for the quantities studied here. Unless otherwise specified, the 
figures refer to an Ecm = ^/s^ = 200 GeV. As we will show at the end of the section, the 
qualitative features do not depend critically on this choice. Furthermore, figures relevant 
for LEP 2 energies can be found in the proceedings of the LEP 2 workshop ^ 



so 



it makes sense to complement here with the higher-energy behaviour relevant for future 
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linear colliders. 



Figure ^ shows the total transverse energy per event for each of the six components 
of the 77 cross section. The spike at small J2E± for the VMDxVMD class comes from 
elastic scattering events, e.g. 77 p^p^. Also diffractive events contribute in this region. 
The large-J^ E± tail of the VMDxVMD curve is enhanced by the possibility of multiple 
parton-parton interactions, which is only included for this class currently. Because of 
the larger Pjfmin cut-off, the classes involving anomalous photons typically have larger 
Y^E±, while the smaller po cut-off for the direct processes corresponds to smaller median 
E^. However, note that the 77 — > qq processes only fall off very slowly with p±, in part 
because of the absence of structure functions, in part because of the form of the matrix 
element itself. The direct x direct class therefore wins out at very large Y^E±. 

The results of Fig. ^ are a bit misleading, since the relative importance of the six event 
classes is not visible. The weighted mixture is shown in Fig. also compared with 7p and 
pp events. One observes a steady progression, with (I]-E_l)pp < {J2E±)^p < {J2E±)j^. 
This pattern, of more activity for a 7 than for a p, is seen in essentially all distributions. 
The elastic spike at small J2 E± is less pronounced for 77, owing to three factors: the 
VMDxVMD component is only a part of the 77 cross section, elastic scattering is a 
smaller fraction of the total p^p^ cross section than it is for pp, and kinetic energy in the 
p° — > 7i~^'7i~ decays add to the total transverse energy. 

The Ej_ flow as a function of pseudorapidity, dEj_/dri, is given in Fig. It illustrates 
how 7p interpolates between pp and 77: around the direction of the incoming photon, the 
7P events look like the 77 ones, while they look more like pp ones in the opposite direction, 
with an intermediate behaviour in the central region. Since elastic and single diffractive 
events are likely to be removed from 'minimum-bias' data samples. Fig. 15 shows the 
behaviour without those two event classes. The quantitative 77/7P/PP differences then 
are shghtly reduced, but qualitatively remain unchanged. Also in subsequent figures these 
events have been removed, and the same comment can be made there. 

The charged-multiplicity distributions follow essentially the same pattern as shown 
for the J2E± ones in Figs. I^^IS, and are therefore not included here. There is one 
noteworthy exception, however: the direct x direct component does not have a tail out to 
large multiplicities. That is, even if the process 77 qq can generate large p± values, 
the absence of any beam jets keeps the multiplicity down. 

The transverse momentum spectrum of charged particles is shown in Fig. |16|. The 
larger high-p^ tail of the 77 processes is one of the simplest observables to experimentally 
establish differences between pp, 7p and 77. Of course, the cause of the differences is to 
be sought in the higher jet rates associated with photon interactions. The jet spectra are 
compared in Fig. Il^, using a simple cone algorithm where a minimum E± of 5 GeV is 



required inside a cone of AR = y (Arj)'^ + (A^)^ < 1. Already for an -E^jet of 5 GeV there 
are more than three times as many jets in 77 as in pp, and this ratio then increases with 
increasing E±jet. The spectacular differences in the jet rate at large p± are highlighted in 
Fig. |18|. They mainly come about because the direct component involves the full energy 
of the incoming photon. The pseudorapidity distribution of jets is shown in Fig. |19|. As 
in the inclusive dE±/dri distributions, the difference in behaviour between the 7 and p 
hemispheres is readily visible. 

To illustrate the energy dependence of these distributions. Fig. ^ gives the dE±/dr] 
flow for cm. energies of 50 GeV. This can be compared with the result for 200 GeV 
in Fig. |15[ Qualitatively, the same pattern is seen at both energies, although relative 
differences tend to be somewhat reduced at larger energies. This is also true for other 
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observables, such as jet rates. One reason is that the possibihty of multiple parton- 
parton interactions in the VMD component pushes up the activity in those events at 
larger energies, and thus brings them closer to the anomalous class. The importance of 
the direct class, on the other hand, is reduced at large energies. Further, at large energies, 
jet production is dominantly initiated by small-x incoming partons, where the VMD and 
anomalous distributions are more similar than at large x (although still different). 

5 Summary and Outlook 

In this paper we have shown that our model for 7p events [§] can be consistently gen- 
eralized to 77 events. That is, essentially all free parameters are fixed by (low-energy) 
7P phenomenology. Since we start out with a more detailed subdivision of the 7p total 
cross section than has conventionally been done in the past, our 77 model also contains a 
richer spectrum of possible processes. We distinguish six main event classes, but most of 
these contain further subdivisions. The aim is that this approach will allow predictions 
for a broader range of observables than is addressed in conventional models. For instance, 
although not discussed in detail here, our approach does correlate the hard-jet physics in 
the central rapidity region with the structure of the beam remnants. 

This does not mean that all results are complicated. We have shown that the simple 
Regge-theory expression (t^J^{s) ~ 211 [nb] comes close to what is obtained 

in our full analysis, and have a fair understanding where differences come from. We 
therefore expect this expression to be good to better than 10% from a few GeV onwards, 
at least to the top 77 energies that could be addressed with the next generation of linear 
e"'"e~ colliders. Also global event properties show a very simple pattern, with more activity 
(transverse energy, multiplicity, jets, . . . ) in 7p events than in pp ones, and still more in 
77 ones. This is perhaps contrary to the naive image of a 'clean' point-like photon. The 
7P events show their intermediate status by having a photon (proton) hemisphere that 
looks much like 77 (pp) events, with a smooth interpolation in the middle. 

In our current model the perturbative approach to the description of the direct and 
anomalous components is pushed to its extreme. In this sense it is a useful study. However, 
we also see that it has its limitations: at high energies a purely perturbative treatment 
of the direct and anomalous components is no longer possible and unitarity corrections 
have to be taken into account, possibly through an eikonalized treatment of the two 
components along the lines indicated above (i.e. treating the direct event class as the 
soft component of the anomalous one). The goal is a consistent treatment covering the 
whole {k±i, k±^2,P±) volume of 77 events in a consistent fashion, with smooth transitions 
between the various regions. Unfortunately this is a not so trivial task, and anyway must 
be based on input from the simpler model above. The current model therefore is a useful 
step towards an improved understanding of the photon and its interactions. 

Also many other aspects need to be studied. Disagreements between the HERA data 
and our model can be found for the profile of beam jets, the structure of underlying events, 
the event topology composition and so on, indicating the need for further refinements |@] . 
The transition from a real 7 to a virtual 7* is still not well understood. Production in 
diffractive systems currently attract much attention. Further issues could be mentioned, 
but the conclusion must be that much work remains before we can claim to have a complete 
overview of the physics involved in 7p and 77 events, let alone understand all the details. 
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Table 1: Coefficients of the total and partial Vp and V1V2 cross sections, according to the 
formulae given in the text, eqs. (0) and (|T9|) . The u is not shown separately, since it is 
assumed to have the same behaviour as the p°. 
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Figure 1: Contributions to hard 7p interactions: a) VMD, b) direct, and c) anomalous. 
Only the basic graphs are illustrated; additional partonic activity is allowed in all three 
processes. The presence of spectator jets has been indicated by dashed lines, while full 
lines show partons that (may) give rise to high-p_L jets. 



Figure 2: The six classes contributing to hard 77 interactions: a) VMD x VMD, 
b) VMD X direct, c) VMD x anomalous, d) direct x direct, e) direct x anomalous, and 
f) anomalous X anomalous. Notation as in Fig. |T[ 
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Figure 3: A generic Feynman diagram (in the leading-log approximation) for 77 interac- 
tions and its decomposition into six components. 



Figure 4: The total 77 cross section. Full curve: the parameterization of eq. ([T0|). Dashed 
curves: range obtained by varying Y^"' as described in the text. Dashed-dotted curve: the 
critical-pomeron parameterization [|T^]. Data points: open triangles PLUTO 1984, filled 
triangles PLUTO 1986, squares TPC/27 1985, spades TPC/27 1991, circles MD-1 1991, 



full square CELLO 1991 P 
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Figure 5: The total VMDxVMD cross section, full curve, and its subdivision by vector- 
meson combination. The components are separated by dashed curves, from bottom to 
top: p°a;, p^J/ip, u!UJ, cucp, ujJ/ip, (jxp, (fiJ/ip, and J/ipJ/ip. Some of the latter 

components are too small to be resolved in the figure. 



Figure 6: The total VMDxVMD cross section, full curve, and its subdivision by event 
topology. The components are separated by dashed curves, from bottom to top: elastic, 
single diffractive (spht for the two sides by the dotted curve), double diffractive, and 
non-diffractive (including jet events unitarized). 
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Figure 7: Components of the 7p cross section: lower dashed curve the VMD contribution, 
middle dashed VMD+direct and upper dashed VMD+direct+anomalous, as obtained by 
integration with ko = 0.6 GeV and the Pj^^(s) of eq. (pll). By comparison, full curve is 
the parameterization of eq. (H). 



Figure 8: (a) Schematic graph for a hard 7p process, illustrating the concept of two 
different scales, (b) The allowed phase space for this process, with the subdivision into 
event classes. 
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Figure 9: Comparison of 77 partial cross sections. Full curves: the simple factoriza- 
tion ansatz of eq. (|2^) . Dashed curves: by integration of jet cross sections (except 
VMDxVMD, where eq. is used). 
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Figure 10: The total 77 cross section. Full curve: the parameterization of eq. ([To|) . 
Dashed curve: result from sum of integrations of the six components. 



Figure 11: Elastic cross sections: p°p° full curve, p^uj dotted, dashed and p^J/ip dash- 
dotted. The latter three cross sections include both mirror-symmetric configurations, 
and have additionally been scaled up by factors 10, 10 and 1000, respectively, for better 
visibility. 
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Figure 12: The total transverse energy per event, separately normalized for each of the 
six event classes. Top frame: VMDxVMD: full histogram; VMDxdirect: dashed one; 
and VMDx anomalous: dash-dotted one. Bottom frame: direct x direct: full histogram; 
direct X anomalous: dashed one; and anomalous x anomalous: dash-dotted one. 
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Figure 13: The total transverse energy per event for different beams: 77: full histogram; 
7p: dashed one: and pp: dash-dotted one. 



Figure 14: Transverse energy flow as a function of pseudorapidity for different beams. 



Notation as in Fig. 13 
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Figure 15: Transverse energy flow as a function of pseudorapidity for different beams, as 
in Fig. |l^ except tliat elastic and single diffractive events have been removed. Notation 
as in Fig. pTS. 



Figure 16: Charged particle inclusive p_L spectra for different beams. Notation as in 

Fig. ra. 
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Figure 17: Rate of reconstructed jets as a function of the transverse jet energy for different 
beams. Notation as in Fig. 



Figure 18: Parton-level jet p± distributions for different beams. The factor 1/2 compen- 



sates for there being 2 jets per event. Notation as in Fig. p. 
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Figure 19: Pseudorapidity distribution of reconstructed jets for different beams. Notation 
as in Fig. 13. 



Figure 20: Transverse energy flow for Ecm = 50 GeV as a function of pseudorapidity for 



different beams, cf. Fig. |T^. Notation as in Fig. [13 
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